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Abstract 

We study inflation with anisotropic hair induced by form fields. In four dimensions, the rele¬ 
vant form fields are gauge (one-form) fields and two-form helds. Assuming the exponential form 
of potential and gauge kinetic functions, we find new exact power-law solutions endowed with 
anisotropic hair. We also explore the phase space of anisotropic inflation and find fixed points cor¬ 
responding to the exact power-law solutions. Moreover, we perform the stability analysis around 
the fixed points to reveal the structure of the phase space. It turns out that one of the fixed points 
becomes an attractor and others (if any) are saddle points. In particular, the one corresponding to 
anisotropic inflation becomes an attractor when it exists. We also argue that various anisotropic 
inflation models can be designed by choosing coupling constants. 
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I. INTRODUCTION 


Scalar fields play a central role in inflationary cosmology. This is because a scalar field 
called an inflaton can mimic a cosmological constant and consequently drive a quasi-de Sit¬ 
ter inflation. On the other hand, other helds such as gauge helds have been overlooked 
in the study of inflation. This is because the cosmic no-hair conjecture insists that other 
matter helds satisfying the dominant and strong energy conditions will be rapidly diluted 
in the presence of the cosmological constant. Indeed, in the case of homogeneous universe. 


fl. 


there exists the cosmic no-hair theorem iJl. Since the inhaton pretends to be a cosmological 
constant, it is natural to expect that this theorem holds even in an inhationary universe. 
Hence, it has been believed that gauge helds are irrelevant during inhation. Nonetheless, 
there have been several attempts for hnding a counter example to the cosmic no-hair con¬ 


jecture by relaxing the energy conditions of matter or by modifying gravity 


S- 


It has been 
. However, recently, a 
Indeed, there exists 


shown that almost all are unstable or need extreme hne tunings 
clear counter example to the cosmic no-hair conjecture was found 
anisotropic inhation with a gauge held hair. The point was not to relax the energy con¬ 
ditions but to notice the slight diherence between de Sitter and quasi-de Sitter spacetime. 


Since the discovery of the anisotropic inhation, there have appeared many phenomenological 
applications jb] in conjunction with the statistical anisotropy j^. 

Because of the phenomenological importance, it is worth clarifying to what extent this 
anisotropic inhation can be extended. In other words, we need to design anisotropic inhation 
as broad as possible. To this aim, we should note that the hrst counter example to the cosmic 
no-hair conjecture consists of a scalar held and a gauge held. From the point of the cosmic 
democracy, however, it is natural to ask if other helds can survive during inhation. Here, 
we do not consider models with higher order derivatives. Rather, we envisage an ehective 
action dimensionally reduced from string theory. As such a model, we consider the action 


S= d‘^x\/—g 


Ml 


R - - n>l>) - E jfVMF, 


p=2 


( 1 ) 


where Mp represents the reduced plank mass, g is the determinant of the metric, R is the 
Ricci scalar, 0 can be regarded as a zero-form held, V (0) is a potential function, and Fp 
are held strength of {p — l)-form helds coupled to an inhaton with gauge kinetic functions 
fp{4>)- In four-dimensions, there are three possibilities, namely, gauge (one-form) helds, two- 
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form fields, and three-form fields. However, three-form helds are equivalent to the potential 
function. In fact, the field strength of three-form is four-form which must be proportional 
to the volume form. Therefore, we do not need to consider this case. One may also think 
two-form fields are equivalent to scalar fields through the Hodge duality = *dxjj with a 
scalar field tp. However, in the cosmological setup, we usually consider homogeneous field 
-^012 (^) = which depends only on a cosmic time. From the view point of a scalar 

held "0, this is an inhomogeneous conhguration ip oc. x^, although this is compatible with 
homogeneity of the universe. Actually, a more complicated version of this possibility is 
known as solid inhation j^. From the point of designing anisotropic inhation, however, it 
is convenient to regard two-form helds as fundamental helds 1^. In this general setup, we 
shall design anisotropic inhation models which exhibit interesting features. 

In this paper, we study inhationary universe in the presence of form helds. In particular, 
we take exponential type potential and gauge kinetic functions. In the case of the gauge 
held with this setup, we know there exists exact power-law solutions a, lUl, ll2|. Hence, we 
expect new exact solutions can be found in the presence of the two-form held. Thus, the 
aim of this work is two-fold. One is to hnd exact power-law solutions in the presence of both 
gauge and two-form helds. The other is to reveal the structure of phase space of anisotropic 
inhation models. 

Since the gauge held and the two-form held give rise to opposite anisotropy, inhation 
becomes isotropic for special coupling constants. While, in general, there arises anisotropy 
of the order of slow roll parameter. Remarkably, co-existence of the gauge held and the 
two-form held gives rise to a complicated phase space structure. Thus, our models provide 
a variety of concrete framework incorporating the statistical anisotropy testable by precise 
observations. 

The organization of the paper is as follows. In section H, we seek power-law solutions 
and obtain exact power-law inhationary solutions with form helds which exhibit anisotropic 
expansion. In section HI, we perform dynamical system analysis and hnd hxed points in 
phase space. It turns out that there are four kind of hxed points corresponding to power-law 
solutions found in section H. We clarify when those hxed points appear. In section IV, we 
examine the linear stability of the hxed points. In section V, we classify the general solutions 
in terms of hxed points. It turns out that only one hxed point becomes stable and the others 
are saddle points if any. The attractor we obtain shows various anisotropy caused by the 
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gauge field or the two-form field. When both gauge held and two-form held survive, we hud 
the convergence to the attractor is slow. We argue how anisotropic inhation can be designed 
by choosing coupling constants. The hnal section is devoted to the conclusion. 


II. EXACT POWER-LAW SOLUTIONS 


In this section, we will seek exact power-law solutions driven by an inhaton 0 in the 
presence of a gauge held and a two-form held For this purpose, we will concentrate 
on the exponential type potential 




and gauge kinetic functions 


/a(0) = , /b(0) = Jb e 


4 > 


( 2 ) 


( 3 ) 


where Vq, /a, Jb are constants. Here, coupling constants X, pa, Pb characterize inhation 
models. Now, we consider the following action 


S= d!^xyj—g 


M, 


fR - :^(a^4>)(d'‘4>) - V(4>) - -Ja\4>)F,,F>‘- - 


( 4 ) 


dyA^ and 


where the held strength tensors and are dehned by = d^ji_Ay 

Hfj,ux = d^Byx + dvB\^ + dxB^y , respectively. It is well known that there exists an isotropic 
power-law inhationary solution for the exponential type potential ([2]). In the presence of a 
gauge held, an anisotropic power-law inhationary solution has been found j^. We shall hnd 
new exact solutions in the presence of additional helds. 

Now, we are in a position to discuss cosmological homogeneous solutions. First of all, 
without loss of generality, we can take the gauge held A^ to be A^ = (0, 0, 0). Given 

this conhguration, in principle, the two-form held can take an arbitrary polarization. How¬ 
ever, in the case of gauge helds, it has been shown that the helds tend to be orthogonal to 
make the expansion as isotropic as possible [h]. Thus, we assume the two-form held B^y 
has an orthogonal conhguration A dx'' = VB{t)dy A dz to the gauge held. Since 

these ansatzes keep the rotational symmetry in the y-z plane, we can take the metric in the 
axially symmetric form 


ds^ = -de + + e^-^^\dy^ + dz^)] , 


( 5 ) 
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where t is the cosmic time. Here, a describes the averaged expansion of the universe and 
a represents anisotropy. From the action (|1]), we can derive the equations of motion. It is 
easy to solve the equations of the gauge field and the two-form held as 


VA=PAffe 


VB =PBfB^^ 2psMpg«+4a ^ 


where Pa,Pb are constants of integration. Substituting these solutions into other equations, 
we can deduce the hamiltonian constraint 


• 2 • 2 I 

a = a + 


1 


3M| , 




.( 7 ) 


the rest of Einstein equations 


1 




1 


2 x-2„-2pA-sT-^-4a-4(T 


1 


a = 




p g 

1 

i 

P 


+ M^Pb/b ^ 


,2 J!-2-‘^PBT$--2a+4o- 


, ( 8 ) 


^Pb/b ^ 


2 ^-2„-2psji^g-2a+4f7 


and the equation of the scalar held 

_L f-2-2pA^-4a-4^. PB 2 r-2 -‘2pB -2a+4a 


= -3d0 - —Hoe + —PaIa e 


M 




/ip g 


To seek power-law solutions, we put the ansatz 
a = C log Mpt , a = T] log Mpt , 


+ jrPBfi e 


/jp g 


Mr, 


= ^ log Mpt + (j)o , 


( 9 ) 


( 10 ) 


( 11 ) 


where (, p, (po are constants. Substituting the above ansatz into the hamiltonian constraint 
o. we obtain 




,2 f-2 -2pA_(po 


6M2 


PaIa e 


(Mpt) 


- 2 pA?- 4 C- 4 p 


+ 


6M2 


2 f-2 -2 pb4'0 


PbIb e 


(Mpt) 


- 2 psC- 2 C+ 4 p 


Taking a look at the power of the time t, we have the relations 

= -2 , 

Pa^ + 2( + 2r] = 1 , 

PB^ + C - ‘2v = ^ ■ 


.( 12 ) 

(13) 

(14) 

(15) 


Here, we assumed form helds are non-trivial. In the case that the held vanishes, the corre¬ 
sponding relation is not necessary. Let us dehne the following variables 


u = —Voe^^° 
M^ 


“ = ^ ■ (16) 
-/Hip 
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Thus, from (fT^ . we have the relation 


0 =-C^+ . (17) 

DODD 

Similarly, from Eqs.(|H]), (E]), and (ITOl) . we can deduce the following relations 


c -3C^ + u + ^w + lz , 

D 6 

(18) 

V — SCp H — w -2) , 

' ^ ' 3 3 

(19) 

^ - 3C^ - uX + WPA + ZPB ■ 

(20) 


A. Isotropic power-law inflation 

First, we start with the well known power-law solution. We consider the case that both 
of the gauge field and the 2-form field vanish, Pa = Pb = 0. Then, Eqs. flTTl) and flTHl) are not 
necessary. Thus, from ffT^ . flTTD ^ fEU]) . we easily obtain the isotropic power-law solution: 

5 = C = ^, -) = 0, a = w = z = 0. (21) 

Notice that we need the condition < 6 for the existence of the solution. For sufficiently 
small A, this describes inflationary universe. 


B. Anisotropic power-law inflation with a gauge field 


Second, we consider the case pa ^ 0 ,pb = 0- 
exact anisotropic power-law inflationary solution 


In this case, it is known there exists an 
j^. The resultant solution is given by 


2 

“a ’ 

(22) 

A^ + 8paX + 12p^ + 8 

(23) 

6A(A + '^‘Pa) 

A^ + ^paX — 4 

3A(A + ^pa) 

(‘^Pa + + 2)(—A^ + 4p^A + 12p^ + 8) 

(24) 

(26) 

2A2(A + 2paY 

(A^ + 2paX — 4)(—A^ + 4:PaX + 12p^ + 8) 

(26) 

2A2(A + 2p^)2 

0 . 

(27) 


6 







From (l26l) . we see the following condition is required for the existence of this power-law 
solution 


2 lx 


(28) 


In inflationary universe A 1, this implies S> 1. 

In order to see the consistency with the cosmic no-hair theorem, we define the averaged 
slow roll parameter 

_ H 6A(A -|- Pa) 

A^ -|- 8p^A -|- 12p^ -|- 8 

The relation between the averaged slow roll parameter and the anisotropy is given by 



E_&_ 2(A2 + 2paA-4) _ 

^ ^ d “ A2 + 8paX + I2p\ + 8 “ 3^ 


where we have defined 

_ A^ -|- 2p^A — 4 

A^ + PaA 

It should be stressed that the anisotropy vanishes in the limit e —)■ 0, which is consistent 
with the cosmic no-hair theorem. 



C. Anisotropic power-law inflation with a two-form field 

Next, we consider the case = 0 7 ^ 0. Then, Eqs. flT^ . flT^ and ffT7|) ^ fl20l) can be 

written as 


Ea.fl^ gives 


From Eq. fl36|) . we have 


Ae = -2, 

(32) 

to 

•Crv 

+ 

1 

to 

03 

(33) 

0 = -c^ + + ^u + ^z , 

0 0 0 

(34) 

0 = C ~ 3^^ + u + -z , 

(35) 

0 = p-3Cp-^z , 

(36) 

0 = ^ - 3C^ -u\ + zpB . 

(37) 


(38) 

^ = -3p(3C - 1) . 

(39) 
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Eqs. fl35|) and (13^ lead to 


u= (3C-1)(C + ^) • 

Using Eqs. fl3^ and fliOj) in Eq. fl37|) . we get 


(40) 


(3C “ i)[ ^ + -^(C + v) + ^PbV ] — 0 • 


(41) 


When C = M and 2 ; become zero. This is not the solution which we desire. Hence, the 
solution should be 


4 + {2pB + \){3pB + A) 
3A(A + pb) 


(42) 


Using this result in Eq. fl33l) . we obtain 


-A2 - 2pbX + 2 
3A(A + Pb) 


(43) 


Substituting the relations fl42l) and fl43l) into Eqs. fl3^ . fl40l) . we reach the solutions 


(6p^ + 4 pbA -|- 4)(2p^ -|- pbA -|- 2) 

(44) 

A2(A + pb)2 

(A^ -|- 2pbA — 2)(6p^ -|- 4 pbA -|- 4) 

(45) 

A2(A + pb)2 


It should be noted that Eq. fl34D is automatically satished. Positivity of Eq. dTSj) requires the 
following condition 

Pb > - - -A . (46) 

If inflation occurs, A <C 1, the above condition implies p^ ^ 1. 

For this solution, the averaged slow roll parameter is given by 


_ 3A(A + Pb) 

A^ + 5pbX + 6p^ + 4 

and the anisotropy is characterized by 

S_ -A2-2pbA + 2 __2^^ 

H A^ + 3pbX + 6p^ + 4 3 

where we dehned 

^ A^ + 2pbA — 2 
2A(A + Pb) 

Again, this is consistent with the cosmic no-hair theorem. 
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(47) 


(48) 


(49) 









D. Anisotropic power-law inflation with a gauge fleld and a two-form fleld 


Finally, we consider the hybrid case that both of the gange fleld and the 2-form fleld are 
non-trivial, pa ^ 0 ^ 0. From Eqs. flTSD ina, we have 



^ ^ + Pa + Pb 

^ ~ 3 A ’ 

_ 1 —A + 2{pa — ‘^Pb) 

^ ~ 6 A ■ 

Snbstitnting these results into Eqs. (fT8|) ^ (l20|) . we obtain 

^ 1 X{2pA + ^Pb) + ‘2{2p\ + -|- PaPb) + 4 

“^2 V 

— 4 + \{2pA — Pb) + 2pb{pa — 2pb) 

VI- ^2 • 

13A^ - 8 + 4A(p^ + Pb) - iPA{PA - 2pb) 

^ ^ 2 A2 


(50) 

(51) 

(52) 


(53) 

(54) 
(56) 


The hamiltonian constraint flT7|) is automatically satisfied. We should mention that the 
inflationary universe is isotropic in the special case, A = 2(p^ — 2pb). In this case, the 
energy density of the gauge fleld is equal to that of the two-form fleld, namely w = z. 
Because of positivity of w and z, Eqs. flSTj) and fl5^ yield the following inequalities 


A^ — 4 -|- \{2pA — Pb) + 2pb{pa — 2pb) > 0 , (56) 

3A^ - 8 -h 4:X{pa + Pb) - ^Pa{pa - 2pb) > 0 . (57) 


In FIGIH we plotted a parameter region for which the hybrid type solution exists. We And 
that as A becomes small, the region becomes narrow and approaches pA = 2pB line. 

Now, the relation between the averaged slow roll parameter and the anisotropy is given 


by 


^ ^ 1 -X + 2{pa-2pb) ^ 

H 4 X + Pa + Pb Q ^ 


(58) 


where 

^ _ 3_A_ ^ _ —A -|- 2{pa — 2pb) 

2 X -\- Pa -\- Pb A 

The consistency with the cosmic no-hair theorem is apparent. 


(59) 
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FIG. 1. We plotted —A+\{2pA — Pb) + ‘^Pb{pa — ‘^Pb) = 0 as the blue curved surface above which 
is an allowed region. The red curved surface 3A^ — 8 + 4A(/9a + pb) — 4:Pa{pa — ‘^Pb) = 0 describes 
the upper bound of an allowed parameter region. The hybrid solutions exist in the parameter 
region between the red and the blue curved surface. 

E. Anisotropic decelerating universe 

Note that we also have exact solutions which represent anisotropic decelerating universe. 
The rate of acceleration is related to the slow roll parameter as 

- = H\1 - e) , (60) 

a 

where a is scale factor and H = a/a . If e > 1, the universe has decelerating expansion. 
Moreover, the relation of the slow roll parameter and the parameter of the equation of state 
w is given by 

w=‘^e-l . (61) 

For example, from Eq.(|29]), we see e can take the value close to 4 for A slightly above 4. The 
corresponding equation of state is tc = 5/3. In the case of two-form dominated solutions, we 
have e = 3 as an upper limit. Thus, it describes anisotropic universe with w = 1. For hybrid 
solutions, we have tc = 0 as an upper limit. This looks like an anisotropic matter dominant 
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universe. So, we can describe the anisotropic decelerating universe with —1/3 < tc < 3 by 
the power-law solutions. 


III. FIXED POINTS IN PHASE SPACE 


In this section, we hnd hxed points in phase space. It turns out that there are four kinds 
of hxed points corresponding to the power-law solutions we saw in the previous section. 

Let us use the number of e-foldings a as a time coordinate. It is convenient to dehne new 
variables 


a 10 

_ _ Y = _— 

a ’ Mr, a ’ 


= jrjPAfA e 
IVlpa 


^-2a-2a 


Zb = 


MpO. 


(62) 

(63) 


With these variables, the equations of motion (|8]) dH can be cast into an autonomous 
system. Using Eqs.dH]) ca, we can calculate the differential of (I6^ ^ (l63|) with respect to 
a as 


^ = A-| 3(A'" - 1) + iyq + Iz^A + 1) + lz|(A - 2) , 


— = Za[ 3(A2 - 1) + M - PaY +1-2X + lzl+ Iz" 

da 2 3 6 

^ = ZBl3(A2-l) + iy2-pBK + 2 + 2A + iz| + lz5 


6 

1], 


(64) 


— (h' + A)[ 3(A^ — 1) + ] + {-Z + Pa + t;)Za + + Pb + t;)Zb > (65) 

dZA 


( 66 ) 

(67) 


The hamiltonian constraint ([7]) is rewritten as 


1 Voe 


\ (P 

At7“ 


M2 d2 


= 3(A= - 1 ) + iy= + Izl + iz| 


Because of the positivity of the potential of inhaton held, the above Eq. 
condition 

0 > 3(a 2 -1) + + izj + iz|. 


( 68 ) 

gives the 

(69) 


We must seek hxed points which satisfy the above constraint. To obtain hxed points, we 
solve Eqs. fl64|) ^ fl671) by setting the right-hand side to zero. The solutions can be classihed 
by Za and Z^- 
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A. Isotropic fixed point: Za = Zb = 0 


First, we look for fixed points in the absence of form fields. From Eg. (16411 . we get 


x [ 3(x'^ -1) + ] = 0 . 


( 70 ) 


Here, 3(A^ ~ 1) + = 0 is against the constraint (16^ . so we have 


A = 0 . 

Thus, from Eg. (16^ . we obtain 

r = -A . 


(71) 


(72) 


This represents isotropic inflation and corresponds to the isotropic power-law solution which 
we saw in the previous section. 


B. Gauge fixed point: Za / 0, Zb = 0 


Next, we take into account the gauge field 



Taking the same procedure as the isotropic 


case, we obtain 


A 

Y 

y2 


2(A^ -|- ‘IpA^ — 4) 

A^ -l- 8pA^ + ^‘^Pa ® 

-12(A + 2p^) 

A^ -|- 8pA^ + -|- 8 

18(A^ -|- 2paX — 4)(—A^ -|- 4pyiA -I- 12p^ -I- 8) 

(A2 + 8pa\ + I2p\ + 8)2 


(73) 

(74) 

(75) 


Note that from fl75|) . this fixed point reguires the condition pA> \ — for the solution to 
exist. And this condition is the same as the condition fl28|) . Also this fixed point corresponds 


to the power-law solution with the gauge field. 


C. Two-form fixed point: 7 ^ 0, Za = 0 


Here, we consider the two-form field. From Egs. (164)1 and (166)1 . we have 


Eg. (16411 gives 


F = -(A + 3pb)A-A . 

2 _ -6A[ 3(A2 - 1) + ly2 ] 


(76) 


(77) 
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Substituting solutions d76|) . fl77|) into Eq. fl67|) . we obtain 


{X + 1)[ (A^ + + ^Pb 4)Ar + + 2pb\ — 2 ] — 0 . 


Apparently, X = — 1 is against the constraint fIBUD . so we take 

^ _ —A^ — 2pbX + 2 

A^ + 5pbX + Qp% + 4 


( 79 ) 


At this time, we have 


Y = 



-Q{X + Pb) 

X^ + 5pbX + 6p^ + 4 
9(A^ + 2pbX — 2)(6p^ + 4p^A + 4) 
A^ + 5pbX + + 4 


(80) 

(81) 


Note that, from the solution flHT]) . we see pb must satisfy the condition pB > ^ — ^A for the 


existence of this hxed point. And this condition is the same as the condition (146|) . 


D. Hybrid fixed point: Za 7 ^ 0, Zb 7 ^ 0 


Finally, we consider the hybrid case. From Fqs. flOOp and fl67)) . we have 

y 4X + 1 _ 

Pb — Pa 

Fq. fl641l gives 

,, 1^2 ^-2 3X[ 3(X^ - 1) + ] 

2^^X + 1 X + 1 

Using Fqs. fl83ll in fl66l) . we obtain 

Zl = 2X[ 3(X2 - 1) + 1^2 I ^ 2(X + 1)[ -3(X2 + pj,Y-1 + 2X] . 
Substituting the solutions fl82l) . fl8^ and fl84D into Eq. fl65|) . we get 


(82) 


(83) 


(84) 


[ (2Pa ~ + Pb +4)X + 2p^ — 2p^ +1 ] [ (4A + 4pyi + 4 pb)X + A — 2p^ + 4p5 ] — 0 . (85) 


For X = _)_4 , the constraint (16^ is not satished. So the hxed point we seek is 

as follows 


1 —A + 2(p^ — 2pb) 
4 A + Pa + Pb 


X + Pa + Pb 

9 A^ — 4 + A(2pA — Pb) + ‘^Pb{pa ~ ‘^Pb) 

4 (A + Pa + PbY 

9 3A^ — 8 + 4A(pa + Pb) — ^Pa{pa — ‘^Pb) 
8 (A + Pa + Ps)^ 


( 86 ) 

(87) 

( 88 ) 

(89) 
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The conditions stemming from the positivity of the Z\ and are the same as the conditions 
fl56|) and fl57|) . This fixed point corresponds to the hybrid power-law solution. 


IV. STABILITY OF FIXED POINTS 

In the previous section, we found the hxed points in phase space. According to the general 
theory of dynamical system, the structure of phase space can be determined by the linear 
analysis around hxed points. Thus, in this section, we examine the linear stability of the 
hxed points which we found in the previous section. 

Now, we linearize Eqs. fl6T|) ^ fl^ around hxed points as 

= 5A[ 3(X2 - 1) + 6X2 ^ 1^2 ^ 1^2 ] 

+<5y[ XF ] + SZaI ^(X + 1)Z^ ]+5Zb[ ^(X - 2)Zb ] , (90) 

= (5X[6X(y + A) ] 

da 

+5y[ 3 (x2 -1) + ^ ^ ^ 1 ^2 ^ 1 ^2 ] 

+5Za[ 2Zyi(-y + Pa + 2) ] + ‘^Zb{-Y + Pb + ^ ] , (91) 

= <5X[ Z^(6X - 2) ] + 5y [ Za{Y - pa) ] 

+dZAl 3(X2 -1) + iy2 _ + 1 _ 2X + zi + ] 

2 6 

+6Zb[ ^ZaZb ] , (92) 

= <5X[ Zb(6X + 2) ] + 5Y[ZBiY - pb)] + <5Za[ ^Z^Zs ] 

da 3 

+5Zb[ 3(X2 - 1) + y 2 + 2X + ^Zl + ] • (93) 

To examine the stability of these equations, we set 

(5X = , 6Y = P^^5Y' , 5Za = e^^SZ'^ , 6Zb = e^^dZ'j^ . (94) 

Then diherential equations reduce to eigenvalue problem. Thus, what we need is to hnd the 
eigenvalues of the system. In the following, we discuss the stability of hxed points case by 
case. 
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A. Stability of isotropic fixed point 


Let us consider the stability of the isotropic fixed point. Substituting the background 
solution Za = Zb = ^-i fITT]) and fl72D into above linearized equations, we obtain 


d{5X) 

da 

d{5Y) 

da 

djSZA) 

da 

diSZs) 

da 


(Ia^ - 3)« . 

(—+ Pa^ ~ 2 ) 32/1 , 
(-A^ + PbA — 1)32b . 


(95) 

(96) 

(97) 

(98) 


The above equations are already diagonalized. Hence, it is easy to read off the eigenvalues. 
Since < 6, this fixed point is stable in the direction of X and Y. On the other hand, 
Eqs.(l97]) and (1^ tell us that the stability requires pA < \ and ps < \ — |A. Taking 
look at the conditions fl28|) and fH6|) . the isotropic fixed point is stable only in the absence 
of other fixed points. 


B. Stability of gauge fixed point 


Next, we examine the stability of the gauge fixed point 


fl. 


The existence of this fixed 


point is guaranteed if the condition 
mat ions A 1, 3> 1 become 

d{5X) 


is satisfied. The linearized equations with approxi- 


da 

d{5Y) 

da 

djSZA) 

da 


= -3SX , 

= -3(5y + V6(A2 + 2paX - 4:)dZA , 
= ~^/6{X^ + 2pAX-4:)6ZA , 


d{SZB) 3(3A^ + 4paX + dp^A - 4p^ + Sp^ps - 8) 


6Z 


B 


(99) 

( 100 ) 

( 101 ) 

( 102 ) 


da A^ + 8paX + 12p^ + 8 

Because of the condition fl28ll . the square root in the above equations is real. Note that X 
and Zb directions are already diagonalized. The eigenvalues are given by 


u; = -3, -- ± *W3(A2 + 2paX - 4) - 


9 3(3A^ + 4p^A + dp^A - 4p^ + SpaPb - 8) 


A^ + 8paX + 12p^ + 8 
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Apparently, the first three eigenvalues have negative real part. So this hxed point becomes 
stable if 3A^ + + dp^A — Ap\ + SpaPb — 8 < 0. This condition is violated when the 

hybrid hxed point appears as one can see from Eq. (l57j) . 


C. Stability of two-form fixed point 


Here, we investigate the stability of the two-form hxed point. Substituting solutions 
and flHTj) into linearized equations, we obtain 

d{6X) 


da 

d{SY) 

da 


= -35X , 

= -36Y + V6(A2 + 2pbX - 2) 6Zb , 




(104) 

(105) 

(106) 


where we used approximations A 1, p^ S> 1 which is the case of inhation. The reality of 
square root is guaranteed by the condition fH6l) . On the other hand, we exactly have 


djSZA) _ 3(A^ + 2paX - Pb^ + ‘^PaPb - 4p| - 4) ^^ 
da A^ + SpsA + 6p^ + 4 ^ 


(107) 


The eigenvalues are given by 


u = 



±*W3(A2 + 2pbA-2) 


9 3(A^ + 2p^A - pbA + 2paPb - ^Pb ~ 4) pc'! 
4’ A2 + 5 pbA + 6p| + 4 


Again, hrst three eigenvalues have negative real part. So this hxed point becomes stable if 
A^ + 2pa\ — Pb^ + 2paPb — dp^ — 4 < 0. Again, this condition is violated when the hybrid 
hxed point appears as one can see from Eq. fl56|) . 


D. Stability of hybrid fixed point 

Finally, we consider the stability of the hybrid fixed point. The equations which we 
obtained by substituting solutions (186|) ^ (189|1 into Eqs. (l90ll ~ (19^ are complicated. So 
this time, we give the results for several representative cases. First, we fix two of the 
coupling constants as A = 0.1, ps = 30. Then, the conditions (l56|) and flFf|) yield 
allowed values 59.9167 < pA < 60.1168. Thus, we take parameter sets {X,pa,Pb) = 
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(0.1, 59.9168, 30), (0.1, 60.0, 30), (0.1, 60.1167, 30) for which the eigenvalues become 


-1.5 X 10“®, 

-3.0, -1.5T3.H) , 

(109) 

-2.3 X 10"^ 

-3.0, -1.5±3.8f) , 

(110) 

-2.0 X 10"®, 

-3.0, -1.5±4.7f) , 

(111) 


respectively. We have also checked other parameters and found the hybrid fixed point is 
stable. Thus, we can infer the fixed point is always stable although we could not prove it 
analytically. Note that first eigenvalues in (11091) ~ fillip are much less than 1. Hence, the 
convergence of this direction is quite slow. 


V. DYNAMICAL STRUCTURE OF PHASE SPACE 



FIG. 2. In the parameter space {X, pA, Pb), we depicted the blue curved surface satisfying the 
relation + 2p^A — pbX + 2pAPB — ^Pb ~ 4 = 0. Below this surface, the two-form field fixed point 
is stable. 
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FIG. 3. We added the red curved surface 3A^ + + 4/9^A — Ap\ + SpAPB — 8 = 0 to FIGi2j 

Above this red curved surface, the gauge fixed point is stable. Note that part of the red curved 
surface is hidden under the blue surface and the same applies to the blue curved surface. 



FIG. 4. The green curved surface A^ + 2paX — 4 = 0 is added to the FIGlS The gauge fixed point 
exists above this surface. 
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FIG. 5. The pink curved surface + 2pbX — 2 = 0 is added to the FIGlU Below the red surface 
and above the blue surface, there are four fixed points of which only the hybrid fixed point is an 
attractor. Above the red surface and right of the pink surface, the gauge fixed point is an attractor 
and the isotropic and two-form fixed points are saddle points. Left of the pink surface and above 
the green surface, the gauge fixed point is an attractor and the isotropic fixed point is a saddle 
point. Below the blue surface and above the green surface, the two-form fixed point is an attractor 
and the isotropic and the gauge fixed points are saddle points. Below the green surface and right 
of the pink surface, the two-form fixed point is an attractor and the isotropic fixed point is a saddle 
point. Below the green surface and left of the pink surface, only the isotropic fixed point exists 
and is stable. 

Now, we are in a position to discuss the dynamical structure of phase space. At the end 
of the day, we will see that there is only one stable fixed point and the others are saddle 
points if any. 

First, we classify solutions in the space of coupling constants. In FIG 121 we plotted the 
blue curved surface satisfying the relation 

+ ^paX — PbX + 2paPb — ^Pb —4 = 0 (112) 

in the parameter space {X,pa,Pb)- Below this surface, the two-form fixed point is stable. 
Conversely, above this curved surface, the two-form fixed point is unstable. Next, in FIGjSl 
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we added the red curved surface 


3A + + 4p^A — '^Pa 4” ^PaPb — 8 — 0 (113) 

to FIG|2l Above this red curved surface, the gauge hxed point is stable. On the other hand, 
below this red surface, the gauge hxed point is unstable. Note that part of the red curved 
surface is hidden behind the blue surface, and vice versa. In fact, below the red surface fl57|) 
and above the blue surface , there exists the hybrid hxed point which we believe stable. 
In FIGJU the green curved surface A^ + 2p^A — 4 = 0 has been added to FIG|31 Note that 
the gauge hxed point exists above this surface (see Eq. fl^H]) ') where isotropic hxed point is 
unstable. In FIGjSl the pink curved surface A^ + 2p^A — 2 = 0 has been added to FIGJU 
In the right of the pink surface, there exists the two-form hxed point. Here, we obtain the 
intersection curve of the green and the pink surface as 

= ( 114 ) 

= ( 115 ) 

Remarkably, substituting these Eqs. (11141) and (11151) into Eqs. (I112p and (11131) . we see these 
four surfaces intersect exactly at the same curve dehned by Eqs. (lll4p and (I115p . 

We can now summarize the results. Below the red surface and above the blue surface, 
there are four hxed points of which only the hybrid hxed point is an attractor and other 
hxed points are saddle. We see that, above the red surface and the right of the pink surface, 
the gauge hxed point is an attractor and the isotropic and two-form hxed points are saddle 
points. In the left of the pink surface and above the green surface, the gauge hxed point is 
an attractor and the isotropic hxed point is a saddle point. Below the blue surface and the 
above the green surface, the two-form hxed point is an attractor and the isotropic and the 
gauge hxed point are saddle points. Below the green surface and the right of pink surface, 
two-form hxed point is an attractor and isotropic hxed point is a saddle point. Below the 
green surface and the left of pink surface, only isotropic hxed point exists and is stable. 
Thus, we can conclude that only one hxed point is stable and others become saddle points 
if any. 

So far, we have investigated the stability of hxed points. Now, we know one is an at¬ 
tractor and the others (if any) are saddle. The trajectory should be very complicated due 
to the saddle points. Any trajectory should approach to the attractor point. However, 
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fig. 6. The coupling constants are set to be A=0.1, /3yi=60.1, pb=30. The number of e-foldings 
is run from 0 to 10. The initial conditions are taken as {X = 0.00015, Y = 0.01, Z\ = 0.07, Zq = 
0.045) for the red trajectory, {X = 0.0001, Y = 0.01, Za = 0.02, Zb = 0.02) for the blue trajectory, 
and {X = 0.0002, y = 0.01, = 0.015, Zb = 0.01) for the green trajectory. The allow represents 

the direction of eigenvector corresponding to the slow mode. The isotropic, gauge, and hybrid 
fixed points are represented by triangle, square, and circle, respectively. The two-form fixed point 
(X, Z^, Zb)=(- 0.000740, 0, 0.0816) is located outside of the range. 


depending on the initial condition, there are several transient periods when the trajectory 
is staying for some time near the saddle point. Here, we show some trajectories for various 
initial conditions. In FIGJHl each trajectory seems to approach an attractor, however, the 
convergent points are slightly different depending on the initial conditions. This comes from 
the fact that the hrst eigenvalues in Eqs. fllOQp ~ fillip are small. Therefore, the trajec¬ 
tories are still on the way to the hxed point. Ultimately, they should converge to the hxed 
point, but it does not reach the hxed point even at the number of e-foldings 100. Hence, in 
reality, we see various anisotropic inhation. We investigated the direction of the slow mode 
numerically (FIG|7p. In practice, the direction of the slow mode through each hybrid hxed 
point becomes an attractor. In FIGJTJ we hxed pb and changed pa in the region where the 
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/.g direction 


X direction 


FIG. 7. The direction of slow mode from each hybrid fixed points. Here, ps = 30, A = 0.1 are 
fixed. Then, we change pA in the region where hybrid fixed point is stable. When pA takes the 
lower limit, slow mode is parallel to the Za direction. This is represented by the red vector. When 
we increase pA-, the eigenvector varies continuously. The vector at pa = ‘^Pb is represented by the 
green vector. Finally, when pA reaches the upper limit, the vector becomes parallel to Zb direction. 
This is represented by the blue vector. 


hybrid fixed point is stable. We also fixed pa and changed pB, and we got a similar graph, 
that is , the direction of the slow mode is determined by the relative ratio between pA and 
2pB- So FI GJ7| represents nniversal behavior of the slow mode. Other interesting trajecto¬ 
ries are plotted in FIG|8l In this case, at hrst, Za and Zb fall to zero, which means the 
nniverse becomes isotropic. Next, each trajectory leave the isotropic saddle point and goes 
to another saddle point where the nniverse is anisotropic. Ultimately, they will converge to 
the stable hxed point where the nniverse nndergoes a different anisotropic expansion. As 
the examples in FIGEJand FIGlHltell ns, we can see that the featnre of anisotropic inflation 
is somewhat inflnenced by the initial conditions. This can be nnderstood from the fact that 
the eigenvalnes of nnstable modes in fllOSp and (11081) are small. Moreover, if we choose a 
different set of conpling constants, we can design a more variety of anisotropic inflation. 
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FIG. 8. We set coupling constants A=0.1, /9 a=140.1, pb=70 and run the number of e-foldings 
from 0 to 15. The initial values are chosen as {X = 0,Y = 0.01, Zj^ = 0.9, Zb = 0.9) for the blue 
trajectory, and {X = 0,Y = 0.01, Za = 0.7, Zb = 0.7) for the green trajectory. The isotropic, 
gauge, two-form, and hybrid fixed points are represented by triangle, square, diamond, and circle, 
respectively. 


Having said the initial condition dependence, one may wonder if the models can be 
predictable. From the FigJTl we see the direction of the slow mode depends on the relative 
ratio between pA and 2pB. Once we hxed the model parameters pa and we can determine 
the eigenvector. Moreover, an initial condition selects a point on this eigendirection. Thus, 
we can read off the magnitude of Za (energy density of gauge held) and Zb (energy density 
of two-form held). Note that Za and Zb are correlated because they are on the line. For 
making predictions, these values are sufficient. Of course, since we cannot know the initial 
condition, we can not predict each observable. Nevertheless, we can falsify the model by 
checking the consistency conditions between observables of the model. In this sense, our 
models have predictability. 
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VI. CONCLUSION 


We studied inflationary universe from the point of view of cosmic democracy. More 
precisely, we investigated phase space of anisotropic inflation in the presence of form fields, 
namely, a gauge held and a two-form held. We sought exact power-law solutions when the 
potential of inflaton and the gauge kinetic functions are exponential type. We found new 
exact anisotropic inflationary solutions. Next, we cast the system into an autonomous system 
and found fixed points corresponding to the power-law solutions. We showed that only one 
fixed point could be stable and the others are saddle points if any. The coupling constants 
X, Pa, Pb determined which one becomes stable. The result is summarized in FIGjS] In 
the case of hybrid fixed point, we found the convergence to the attractor is very slow. As 
a consequence, it turns out that the dynamics of the anisotropic inflation depends on the 
initial conditions. We also found trajectories which go around several saddle points before 
reaching the final attractor point. We also demonstrated how the phase space of anisotropic 
inflation is fertile. 


There are several applications to be considered. Our results may have implication for 
generation of primordial magnetic fields 13| . It is worth investigating generation of magnetic 
fields in this context of cosmic democracy. In the present model, an isotropic inflation can 
be realized in spite of the presence of the gauge held and the two-form held because both 
tend to produce opposite anisotropy (This is similar to the multi-vector cases InlQ.). 
It is intriguing to study observational signature of this special case. In general, of course, 
we have anisotropic inflation. It is also interesting to look at observational predictions of 
hybrid models on the cosmic microwave background radiation. Depending on the solutions, 
the predictions would be different. For example, the shape of the anisotropy depends on the 
initial conditions. However, that does not mean the model is not falsiflable. In fact, there 
are several consistency conditions between observables which enable us to get information 
independent on initial conditions. We leave the detailed study to future work. 


From the perspective of designing anisotropic inflation, there are several directions to 
be investigated. It is possible to extend anisotropic inflation to non-abelian models 15| . 


k-inflation models 
tion 


16| . We can also incorporate the parity violation in anisotropic infla- 


ITj . These are complementary to the effective held theory approach 181 
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